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1 INTRODUCTION 

In a cross field device[|lj such as magnetron or cross field 
amplifier electrons move in crossed magnetic and electric 
fields. Due to synchronism between electron drift velocity 
and phase velocity of RF wave, the wave bunches the beam, 
electron spokes are formed and the bunched electrons are 
decelerated by the RF field. Such devices have high effi- 
ciency (up to 90%), high output power and relatively low 
cost. Electrical design of the cross-field devices is diffi- 
cult. The problem is 2D (or 3D) and highly nonlinear It 
has complex geometry and strong space charge effects. Re- 
cently, increased performance of computers and availabil- 
ity of Particle-In-Cell (PIC) codes ^], have made possi- 
ble the design of relatively low efficiency devices such as 
relativistic magnetrons or cross field amplifiers [^. Sim- 
ulation of high efficiency 90%) devices is difficult due 
to the long transient process of starting oscillations. Use 
of PIC codes for design of such devices is not practical. In 
this report we describe a frequency domain method that de- 
veloped for simulating high efficiency cross-field devices. 
In the method, we consider steady-state interaction of par- 
ticles with the modes of RF cavity at dominant frequency. 
Self-consistency of the solution is reached by iterations un- 
til power balance is achieved. 

2 PHYSICAL MODEL 

Cross-field devices consist of a cathode and a surround- 
ing anode. The structure is a cavity with a set of resonant 
eigenmodes. Macroparticles are emitted from the cathode 
and moved by forces of electromagnetic fields. The elec- 
tromagnetic fields are determined by applied external elec- 
tric potential between anode and cathode, oscillating field 
of cavity modes, and space charge fields. We use geom- 
etry with arbitrary piece-wise planar boundaries. In order 
to solve the electrostatic and electrodynamic problems, we 
apply methods that do not require mesh generation. Inter- 
action with magnetic field is determined by uniform mag- 
netic field Hz which is parallel to z-axis and orthogonal 
to the plane of simulation. There are several assumptions 
that we use to simplify the problem. These assumptions are 
based on the working regime of the devices that we want to 
simulate. Devices will have low current density, are non- 
relativistic, and have resonant systems with a relatively low 
density of the cavity modes. Hence, we can neglect mag- 
netic fields due to space charge and cavity modes. We can 
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also use cavity modes with eigen-frequencies close to the 
working frequency. 

2. 1 Basic equations 

We are solving a steady state problem of electron beam 
flow in self-consistent electromagnetic fields. Total fields 
are superposition of static electric E' and magnetic H' 
fields, and "oscillating" electric E{uj) and magnetic -ff (w) 
fields as 

E{t) = #+Ke{i?(w)e^"*}, H{t) = H'+die{H{uj)e^'^*}. 

Here cu is angular frequency, t is time. We separate the elec- 
trodynamic problem into two parts. The first part - electro- 
static potential $ is generated by "external" anode-cathode 
potential and by the static component of the space-charge 
electric fields. The second part - the dynamic electromag- 
netic fields have a harmonic e^'^* time (t) dependence. 

2.2 Static fields 

We find the static electric field from E' — — V$, using the 
Poisson equation: 

S/H = -^, (1) 

where V is the gradient operator, p is volume charge den- 
sity averaged over oscillation period T ~ 2tt/u!. eq is the 

electric permittivity of the vacuum. 

2.3 Oscillating fields 

To solve the second part of the problem, we write the time 
harmonic Maxwell equations as 



V X iJ = -jujpoH, V X H = jujeaE + 



(2) 



Here J^j is electric current density, is the magnetic per- 
meability of vacuum. Oscillating fields inside a cavity are 
expanded in terms of the cavity eigenmodes (E^, Hg) and 
the fast oscillating electric potential ip^ as 



E = Y, AsE, - V(^^, H = Y1 B^^s. 



(3) 



Here s is mode index. As and Bg are the eigenmode am 
plitudes. Using the expansion (^ we get the Poisson equa 
tion for the potential: 

^ ■ Jlu Puj 



V> = 



jLueo 



(4) 



where p^^ is the oscillating space-charge density. Ampli- 
tudes of the electric field expansion are given by 



A, 



jiuj^-uj^s)eoJyEgE*dV 



(5) 



Here ujs is the mode eigen-frequency of the mode , V is the 
cavity volume. 



2.4 Equation of motion 

Equation of motion for an electron in crossed-fields is 



dip 
~dt 



(6) 



where p is the relativistic momentum, is the charge, and 
V is the velocity of the electron. Current density induced by 
the electron motion is J = qevS{r), where r is the position 
vector of the electron, and 6 is the Kronecker delta function. 

3 NUMERICAL METHODS 

We created several separate program modules to simulate a 
cross-field device. First is an RF field solver that calculates 
eigenmodes and eigen-frequencies in the cavity; second is 
the Poisson solver that finds electric fields due to external 
potential, static space charge, and oscillating space charge; 
and third, the tracking module that performs tracking of 
electrons through electromagnetic fields. For simulation, 
we consider an arbitrary, piecewise bounded 2D geometry. 

3. 1 Planar geometry 




Figure 1 : Planar geometry. 

The geometry is cylindrical (uniform in the z-direction) 
as illustrated on Fig. |l]. It consists of planar sidewalls and 
apertures. The geometry in the x, y plane can be described 
by a set of points Zg = (xs, ?/s), where s = 1, 2..., A^'; here 
N' is the total number of sidewalls and apertures. Periodic 
boundary conditions are applied to the apertures. The pe- 
riodic boundary allows us to use only part of the structure 
and significantly reduce simulation time. In the particular 
case shown in Fig. 0the geometry has A^' = 19 sidewalls, 
two apertures (ports) with starting points p ~ 1, 15, and 
the cathode and anode determined bys — 16, 17, 18, 19, 1 
and s = 2,3,...,15 respectively. 

3.2 RF field solver 

The description of the RF solver that is used in this method 
is published in [||]. Here we briefly outline its properties. 
We use the scattering matrix approach to calculate the 
dispersion parameters of the periodic 2D structure, it's res- 
onant frequencies, and the corresponding fields. The fields 
are described by functional expansion. Boundary contour 
mode-matching is applied in a piecewise bounded 2D re- 
gion is applied to obtain the scattering matrix and field am- 
plitudes [0]. The Galerkin method is used for the mode- 
matching procedure. The geometry is divided into regions, 
and electromagnetic fields in each region are expanded in 



series of plane waves or (for low frequencies) Bessel func- 
tions. Scattering matrices from the regions are combined 
using the generalized scattering matrix technique. Reso- 
nant and periodic boundary conditions [^j] are used to ob- 
tain resonant frequencies, dispersion parameters, and cor- 
responding fields. We calculate the electric fields on a polar 
grid (only in the region of field-particle interaction), in or- 
der to speed up calculation of fields for the macroparticle 
tracking. To obtain field at the macroparticle position we 
use 2D spline interpolation. 

3.3 Poisson solver 

We use an efficient method for solving the Poisson equa- 
tion for electric fields in a 2-D, arbitrarily shaped geometry. 
The solution is based on the method of moments. Point- 
matching in a piecewise bounded 2D region is applied to 
obtain the charge density on the boundary. The boundary's 
charge density determines the fields and potentials through- 
out the interior region. We use a complex representation 
of the fields and potentials in the solution [^. We apply 
periodic boundary conditions to simulate the fields in the 
periodic structure. 

Formulation We solve equation (|l|) in 2D. In the 2D 

case it is advantageous to represent the position and field 
vector's (x, y) components by a single complex represen- 
tation. We will work with functions of a complex variable 
z = x + jy. The field strength E can be written in terms of 
the scalar potential $ = ^{z) as 



Eiz) 



dz 



(7) 



Here * represents the complex conjugate. An effective line 
charge q (point charge in 2D geometry) has the complex 
potential $ — (<z/eo) log z. We approximate the charge dis- 
tribution on the boundary of the region as a sum of "step" 
functions. We divide each element(sidewall and aperture) 
of the boundary into Nf, straight pieces or "charged lines" 
with uniform charge density a along the piece. A uniformly 
charged straight wall with beginning and end coordinates 
zi and Z2, respectively, will produce a complex potential at 
the point z^ 



L eo 



■ log(z - Zio)dz, 



(8) 



where L is the contour along the line. Equation ^ is inte- 
grated analytically. 

Field strength of the charged wall We obtain the 
electric field of the charged line by substituting (^ into (^: 



E{z^)eo 



\Zl - Z2\ 
Z\ - Z2 



log 



Zyj Z\ 

Zw - Z2 



(9) 



The value of the function is undefined on the line's contour 
However, for us, the fields inside the region are of interest. 
Therefore, the direction of the field (for positive charge) 



on the line's contour is chosen to be directed inward. Also 
singularities at points zi and Z2 can affect the field's cal- 
culation. Macroparticles with finite dimensions are used to 
avoid this singularity. 

Periodic boundary condition We assume that the po- 
tential and field strength are repeated on the period's aper- 
hires (Fig. [l|). Let z[ e Y( and e Y^. If we shift the 
region to the right so it coincides with the next period, the 
coordinate will be transformed into coordinate ^3. The 
periodic boundary condition becomes 



dn 



dn 



We assume the Dirichlet condition on the sidewalls (except 
for the apertures) as 



$(r') = c(r'), r = r' 



(11) 



Integral equations For periodic boundary conditions 
( p^ and dll] ) surface charge density cr must satisfy the cou- 
pled integral equations 



/plog(z^ - z)a{z)dz = eoCizw), e T' , 
/p \og{z[ — z)a{z)dz — /p log(z2 ~ z)a{z)dz, 



aiog(4-z) 



dz 



I (j{z)dz — TTa{z2), 



zeT,z[eY{, z'^eYi, 



(12) 

in which ^ iog(^™ ^) (jgnotes the normal derivative of 

orip 

log(z^, — z) at the point 2^, assuming z is fixed; T = 

r' + Y( + Y^ 

([To|); and C{^io) is the external potential 



Fn'; coordinates zi and Z2 are the same as in 



Numerical approximation We solve the integral 
equation numerically, by approximating the source densi- 
ties by step-functions |^ . Thus we divide the given bound- 
ary r into iVr intervals and assume that the simple source 
density a has a constant value within each interval. Then 
denoting these constant values by di, i = 1,2, Nr, we 
approximate <i> and E by 

^{zw) = ^cr^eo / ^og{z^ ~ z)dz, and (13) 

,■1 



i=l 



(Tt f f d\og{z^ - z) 



dz,. 



dz, (14) 



where denotes integration over the i-th interval of F. We 
substitute ( |l3l ) and (14) into (12) to obtain numerical ap- 
proximation for periodic solution. The unknowns (in the 
system obtained) are the charge density on the intervals tJi, 
the potential and the electric field on the periodic aperture. 
All coefficients in the system are calculated analytically. 
For practical geometries, the matrix of coefficients is well 



defined and there is no difficulty in solving the system di- 
rectly. For macroparticle tracking, the electric field calcu- 
lated on polar grid and then interpolated at the macroparti- 
cle position (same as for RP fields). 

3.4 Tracking 

We find a macroparticle trajectory by using the 4th order 
Runge-Kutta method for integrating the equation of mo- 
tion (H) in polar coordinates. Then, we integrate the com- 
plex electric field of the cavity modes along the trajectory to 
find coefficients for the cavity's eigenmodes (|^). We mon- 
itor energy conservation in order to verify accuracy of cal- 
culation. For that purpose we use total energy that consists 
of kinetic energy of the macroparticle and integral of static 
(due to external potential and static space charge) and os- 
cillating (due to cavity modes and oscillating space charge) 
electric fields along the trajectory. Initial charge and ve- 
locity V are determined by a space-charge-limited-emission 
model and a relaxation scheme. 

3.5 Algorithm 

We start simulation by calculating dispersion the curve for 
the spatial period of the device (using the RF field solver). 
Then, we calculate electric fields for the eigenmodes. Next, 
(using the Poisson solver) we calculate electric field due to 
external potential. Next, we start iterations using Track- 
ing module to find the macroparticle trajectories, field inte- 
grals along the trajectories, and electric fields due to space 
charge. Next, we update the static and oscillating fields and 
start new iteration. 

4 SUMMARY 

We have written a C-n- computer code that uses meth- 
ods, described above. Accuracy of resonant frequency 
calculation by RF field solver for typical geometries is 
~ 0.1%. We tested performance of Poisson solver and 
Tracking module on diode geometries (without magnetic 
field). We calculated diode current with typical accuracy 
2-3% in comparison with analytical solution. Testing of 
the code on cross- field devices is under way. 
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